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Abstract. We measure the variance in the distribution of off-plane (|6| > 20°) galaxies with rriK < 13.5 from the 
2MASS K-band survey in circles of diameter between 0.344° and 57.2°. The use of a near-infrared survey makes 
negligible the contribution of Galactic extinction to these fluctuations. We calculate these variances within the 
standard ACDM model assuming that the sources are distributed like halos of the corresponding mass, and it 
reproduces qualitatively the galaxy counts variance. Therefore, we test that the counts can be basically explained 
in terms only of the large scale structure. A second result of this paper is a new method to determine the two 
point correlation function obtained by forcing agreement between model and data. This method does not need 
the knowledge of the two-point angular correlation function, allows an estimation of the errors (which are low 
with this method), and can be used even with incomplete surveys. 

Using this method we get f(z = 0,r < 10 /i _1 Mpc) = (29.8±0.3)(r//i~ 1 Mpc)~ 1 ' 79±0 02 , which is the first measure 
of the amplitude of £ in the local Universe for the K-band. It is more or less in agreement with those obtained 
through red optical filters selected samples, but it is larger than the amplitude obtained for blue optical filters 
selected samples. 

Key words. Large-scale structure of Universe — Infrared: galaxies — Galaxies: statistics — Cosmology: theory 



1. Introduction 

The distribution of the number of galaxies or clusters 
of galaxies in a certain volume V can be used to test 
the large-scale structure. The calculation of probabilities 
of different structures in a clustering model are some- 
times amenable to analytical techniques: for instance, 
for voids (Betancort-Rijo 1990), density of Abel clusters 
(Betancort-Rijo 1995), number of galaxies in a randomly 
placed sphere and cluster density profiles (Betancort-Rijo 
& Lopez-Corredoira 1996), etc. The calculation of variance 
of the number of galaxies within the cones corresponding 
to the circles in the sky, or another 2-D figure, (i.e., galaxy 
counts) is another possibility, and even easier to compare 
with the observations since we do not need any informa- 
tion about the redshifts of the individual galaxies. This is 
precisely the purpose of this paper: the calculation from 
a model of the variance of the number of galaxies in the 
different regions of the sky, and to see how well they fit 
real data obtained directly from a sky survey. This will 
allow us to determine some parameters of the two-point 
correlation function independent of any other method. 

One disadvantage of the galaxy counts is that they 
are affected by Galactic extinction, and it is difficult to 
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separate the fluctuations due to this extinction from the 
real fluctuation in the galaxies distribution. However, the 
arrival of near-infrared galaxy surveys such as DeNIS or 
2MASS provides us a way to avoid this problem, since the 
extinction in K-band, the one to be used in this paper, is 
10 times lower than in V-band. In V-band, the dispersion 
of galaxy counts due only to the fluctuations of extinction 
in scales around 1 degree is around 30% [see eq. (@J and 
the parameters given in the paragraph before eq. flf] at 
galactic latitudes \b\ > 20°, while in K-band is only ~ 3%. 
Since the fluctuations we measure due to the large-scale 
structure are 25-80% in the explored scales, it is clear that 
the calculation of the structure parameters are very sen- 
sitive to the exact knowledge of the galactic extinction in 
visible bands; however, in near-infrared, the calculation of 
the parameters is much less sensitive to the extinction and 
we can be sure that the measured fluctuations correspond 
to the real distribution of galaxies. Moreover, in the K- 
band we have an important advantage: the K-correction is 
nearly independent of galaxy type (Mannucci et al. 2001), 
and this will allow us to make the K-corrections without 
knowing the galaxy type of the sources in the survey. 

Maps of galaxy counts have been available for quite a 
long time and have been used for several purposes. Zwicky 
(1957, p. 84), for instance, produced them from the avail- 
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able catalogues in visible bands, with a smooth correction 
of galactic extinction as a function of b, in order to ex- 
plore the intergalactic extinction. Rather than obtaining 
the Galactic or intergalactic extinction, the purpose of this 
paper is to show that the count variances are in qualitative 
good agreement with those obtained assuming that the 
source correlation function is equal to that obtained in the 
standard ACDM model. Cumulative galaxy counts (in op- 
tical) versus the limiting magnitude were used to constrain 
the departure from homogeneity — fractal distribution — at 
large scales (Sandage et al. 1972). A log N(m) cx 0.6m 
points out total homogeneity and this is more or less ob- 
served in near infrared counts too: with 2MASS data for 
low redshift (Schneider et al. 1998, fig. 5) or other surveys 
(e.x., Kiimmel & Wagner 2000). Here, we will not explore 
further the homogeneity at very large scales but the clus- 
tering at small scales for the local Universe (z < 0.4). 
The use of galaxy counts to derive the two-point angu- 
lar correlation function was already used by Porciani & 
Giavalisco (2002); however, we will recover directly (un- 
der some assumptions) the two-point correlation function 
of the sources, and, therefore, the relative biasing with 
respect to the mentioned halos. We may accurately deter- 
mine it, by considering a slightly modified power law (the 
correlation function is a power law in the most relevant 
region) and determining the value of the amplitude and 
the exponent which lead to the best agreement. 

The paper is divided as follows: ^describes ^ e obser- 
vational data and the way to measure the variance in the 
distribution of galaxies; explains how to calculate this 
variance within a model for the large scale structure; and 
2] makes the comparison between data and model pre- 
dictions and derive the parameters of the two-point cor- 
relation function necessary to get an agreement between 
model and data. 

2. Galaxy counts 

The data used in this work have been taken 
from the extended sources of the 2MASS- 
project (Jarrett et al. 2000), All-sky release 
(http:/ /www. ipac.caltech.edu/2mass/releases/docs. html). 
Completeness limit: m,K = 13.5 (Schneider et al. 1998; 
Jarrett et al. 2000; Mailer et al. 2003, §2). Assuming 
an average color of B — K w 4, it is equivalent to an 
optical limit of 17.5 (Schneider et al. 1998), deep enough 
for statistical studies of the large scale structure, the 
local structures are not too predominant in it, although 
high redshift are excluded [the galaxies have redshifts 
z <0.3-0.4 (Cole et al. 2002), and an average (z) w 0.083 
(see below for details)]. We do not analyze samples of 
galaxies below this limit of ttik (for instance, mx < 12.0 
or rriK < 11.0, etc) because this would represent the very 
local Universe rather than the large scale structure. 

In Fig. ^ we see a representation of the fluctuations at 
scale of 3° in the whole sky (average around 150 galaxies 
per area with complete coverage). In the zone of avoid- 
ance, |6| < 20°, there is a clear defect of galaxies due to 



the extinction, small compared to the optical but not neg- 
ligible in near plane regions. The reliability is larger than 
99% in \b\ > 20° (Schneider et al. 1998). 

In order to quantify the fluctuations, we count the 
number of galaxies in each circle of sky with angular ra- 
dius ro. In 3D space, it counts the galaxies within the 
corresponding cone in the line of sight. We select only 
the circular regions in |6| > 20° which were covered in 
more than 90%. The cumulative counts up to magnitude 
13.5 are expressed per unit area (we divide the number of 
galaxies per region by the area, S of the region). We mea- 
sure the counts in randomly placed circular areas instead 
of a regular mesh. Since the number of random circles in 
which we measure the galaxy counts is around 8 times 
larger than the total area divided by the area of the cir- 
cle, the lost of information is very low (a 0.033% of the 
galaxies would not be in any circle if the distribution were 
Poissonian). Once we have these counts for each of the n 
regions containing respectively rii galaxies, we calculate 
the average, 
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14.4 deg , (1) 
and the dispersion with respect to the average, 



cjv(0o) 



\ 




(2) 



This last number gives us information about the amount 
of structure in the 3D cones with angular radius 9q. An 
example of the distribution of counts is given in Fig. [5] 
and compared with a negative binomial distribution with 
the same variance, to which a wide variety of clustering 
processes lead (Betancort-Rijo 2000). 

The fluctuations due to the intrinsic clustering of 
galaxies, <7 s t, will be the total fluctuations minus the other 
independent sources of fluctuation subtracted quadrati- 
cally: 



(SN) a 



N 



Poisson ' 



(3) 



The galactic extinction r.m.s. fluctuations, cr cxt , may 
be evaluated using the Schlegel et al. (1998) maps of ex- 
tinction: in |6| > 20°, the mean extinction in K-band is 
Ak — 0.019 mag, and the fluctuations of this extinc- 
tion are <ja k ~ 0.020 mag (in the scales with 9o around 
1 degree; it changes slightly with the scale). Given that 
log iV = constant + Q.Qvtlk (for a homogeneous distribu- 
tion; see, for instance, Sandage et al. 1972): 



0.6(lnlO)o- Aj< < N > 



(4) 



For the K-band <7jv >> c cxt , so the small uncertainties 
in the knowledge of the statistical properties of the ex- 
tinction are not important. Note however that in visible 
bands <7 ex t will be 10 times larger, comparable to the dis- 
persion due to the large scale structure and, therefore, it is 
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Fig. 1. Aitoff projection in galactic coordinates (up) and equatorial coordinates (down) of the 2MASS/ "All-sky release" 
galaxy counts with ttlk < 13.5; average: 14.4 galaxies/deg 2 . Galaxies were counted in square regions of 3 deg.x3 deg. 
A smoothing of the map was carried out. Black color stands for regions with very few galaxies, and high density is 
represented by white color. Extinction effects are negligible in the fluctuations. Poissonian noise is not negligible but 
it is relatively small compared to the intrinsic fluctuations due the large scale-structure (see text). 



not possible to obtain accurate information about this last 
dispersion without an accurate measure of the extinction 
fluctuations. This is precisely the advantage of infrared 
surveys. 

The Poissonian fluctuations is the most important con- 
tribution, apart from those due to large scale structure, 
and it can be exactly determined: 



<j Poissou ee ((SA0 2 >PoL on = V<N><1/A>. (5) 

Once the fluctuations due to the structure are obtained 
by means of eq. 10 , these can be compared with the pre- 
dictions of a standard model, as it will be carried out in 
the following sections. 



3. Model 

The calculation of the fluctuations ^ = v ^'^ 3t due to 
the clustering of the galaxies is carried out in the following 
way: 

(AO = An [ rH[M K < M KMmit (r)]dr, (6) 

{5N)l=n 2 [ dn f $[M K < M K ,timtt<Ti)] 

x $[M K < M*r ; ; imit (r 2 )]£[|ri - r 2 \, z(n)]dr 2 , (7) 

where r stands for the radial distance in comoving coordi- 
nates, A = 27r(l — cos#o) is the area (in steradians) of the 
circular regions, n is the mean space density of galaxies 
(which is irrelevant in this context since it cancels when we 
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Fig. 2. Distribution of frequencies of cells with N galax- 
ies for circular cells of angular radius 6q = 0.040 rad. 
Binning step: 10. For comparison, we plot a negative bi- 
nomial distribution (Betancort-Rijo 2000) with the same 
average and variance than the data. 

calculate SN/N), Q[Mk < MKAtmitir)] is the cumulative 
normalized luminosity function up to absolute magnitude 
M K ,iimit[r P hys = r/(l + z)] = m K ,Umit + 5 -b\ogr phys - 
CorrK(r p h ys ). The K-correction is CorrK[r p h ys = + 
z)] = -2.955z(r phys ) + 3.321z(r p/l!;s ) 2 for z < 0.40 (be- 
yond z = 0.40 we do not see practically any galaxy), 
for any galaxy type (obtained from the fit of data from 
Mannucci et al. 2001). The redshift is z(r p h ys ) corresponds 
to the velocity v{r p h ys ) — 100r p ^j )S (/i~ 1 Mpc) km/s (lin- 
ear Hubble law taken as an approximation for low red- 
shift galaxies). The luminosity function $ is taken from 
Kochanek et al. (2001). We neglect the evolution of the 
luminosity function, which is mild ( dl ^ z K ~ 0.5, Pozzetti 
et al. 2003) and very small at the mean redshift of our 
galaxies ((z) as 0.083): a shift of ~ 0.04 magnitudes will 
produce a variation of ~ 2% in both (SN) and (N) in the 
same direction; since we calculate the rate ^jf, this small 
variations mostly compensate, and only second order vari- 
ations remain, which are negligible. The errors associated 
with the uncertainty in the luminosity itself do mainly af- 
fect the normalization, which is also irrelevant for the cal- 
culation of SN/N since normalization factor cancels out. 
Finally, the two-point correlation function is taken to be 
(for comoving coordinates): 



£(r,z)=£(^ = 0)(l + z r (3+e+ ^, 



■1, 



£(r, z = 0) = 



r < 0.01 
0.01 < r < 10 



(8) 



(9) 



FT[P(K)], r > 10 



where FT[...] stands for Fourier transform, and r is given 
in units ft, _1 Mpc. The truncation of the two-point corre- 
lation function is set at 10 /i _1 kpc, but the exact value 



does not matter, our results do not depend on this mini- 
mum scale; it is introduced to avoid computer algorithm 
problems at very small scales. More exactly, the change of 
regime should occur at the matching point but, for the in- 
teresting range of parameters, this takes place around 10 
/z -1 Mpc. The power spectrum for the linear part (r >ss 10 
/z -1 Mpc), P(K), is given for a CDM scenario with initial 
Harrinson-Zeldovich power spectra (Bardeen et al. 1986): 



P(K) 
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A\n(l + 2.34qf 



Ky/l + 3.89g + (16. lq) 2 + (5.46g) 3 + (6.71<?) 4 



K 



n mat h (/i/Mpc) : 



(10) 



and we adopt A = 1.984 x 10 4 (er 8 = 1) and Q ma .th = 
0.21. Within small scales (9q < 0.040 rad), the non-linear 
part is predominant, so slight variations of the parameter 
^mat^i or the normalization A will not lead to changes 
in our results. For instance, numerical experiments have 
shown us that a variation of 10% in the amplitude A leads 
to variations of w 2% in SN/N(6 = 0.040 rad) or w 
0.7% in SN/N(6 = 0.010 rad) and much less for lower 6> ; 
since the error bars of 5N/N are of this order, we consider 
negligible these small variations due to the change of the 
parameters in P(K). In scales larger than 6q — 0.040 rad, 
the errors bars are too large compared to the possible 
variations in P(K). In all ranges, even variations up to 
10-20% in P{K) will not affect SN/N too much compared 
to the errors. 

Two questions might rise as to the suitability of eq. (JjJJ : 
1) is the assumption of a power-law for non- linear scales 
appropriate?; 2) could we apply a power-law for all scales 
instead of a ACDM model for the non-linear regime? Both 
questions are answered in other papers but they will also 
be answered by the result of the fit of the counts itself, 
shown in jJJI (see Fig. OJ. As to the first question, the fit 
of the power law with ro = 6.66 h~ 1 Mpc, 7 = 1.79 is 
remarkably good for low 9q (which is nearly independent of 
the linear part of £), a power law in the non- linear regime 
gives a very good fit. The answer to the second question 
is provided by the following consideration: a power-law in 
all scales gives more structure at r > 40 ft, _1 Mpc than the 
power-law in the non-linear regime + ACDM model in the 
linear regime. The first option gives further fluctuations 
at large 9q] the difference is not high enough to reject the 
first option (we have already said that the fit is not very 
sensitive to the parameters in the linear power spectra), 
but the ACDM model in the linear regime (solid line in 
Fig. is considerably better than the power-law at all 
scales (long dashed line in Fig.|3). 

There are also theoretical reasons to use expression iJJJJ 
for £(r). ACDM simulations lead to a correlation function 
for virialized halos with circular velocity > 120 km/s with 
this form and 7 = 1.7 (Primack 2001, Fig. 1). This func- 
tion turns out to be almost exactly equal to that for APM 
galaxies. For K-selected galaxies the correlations do not 
need to be equal to that for halos; there may exist some 
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biasing. However, the biasing (b{r) — f g~~[7j ) ) is ex- 
pected to be a mild function of r. So, it seems plausible 
to use expression (0 with different amplitude and slightly 
different 7 with respect to the halos (or APM galaxies). 

With this in mind, we shall see in next section that 
with the assumed shape in eq. @, there is qualitative 
agreement with the observed variances. We consider this 
result as a confirmation of the assumption that the mea- 
sured variances are due to the large scale structure. This, 
along with the previously given arguments supporting this 
assumption, make us confident as to its correctness. So, we 
consider that large scale structure should explain not only 
the main part of the variances, but the whole of them. 
Forcing this by choosing the appropriate values for ro, 7 
lead us to an alternative method for determining the cor- 
relation function. 

The evolution of the correlation function depends on 
fi mat and Oa through e. In our case, since the average red- 
shift (z) = 0.083 (for itik < 13.5 using the aforementioned 
luminosity function and K-correction) , the exact value of 
e is not so important, and small variations on it will not 
affect significantly the results from the model. The evolu- 
tionary corrections in this small range of redshifts are also 
negligible, specially in K-band (Carlberg et al. 1997). We 
take the value e = —0.1, which comes from the approx- 
imation of £ cx D(z) 2 , for comoving coordinates [which 
holds provided that the shape of £ does not evolve, i.e. 7 
is constant with respect to z (proved by Carlberg et al. 
1997)], where D, the growing factor of the linear density 
fluctuations, is (Heath 1977, Carroll et al. 1992) given by 



D(z) = 



R = 



5n mat (t )(i + z) 

2f[R(z)} 



f(R)dR, 



(11) 



f(R) = 



I +fimat(*o) ( ^ - 1 



n A {R 2 -i) 



-1/2 



This leads to a value of e « —0.1 for r2 mat (io) = 0.3, 
£1a = 0.7. Although, as we have said before, this value 
changes with other cosmological parameters, our outcome 
will be practically independent of these small changes, e 
was measured by Carlberg et al. (1997) using a K-band 
deep survey, and although the accuracy was not very high, 
they obtained a value (e = 0.2 ± 0.5) which is compatible 
with e = -0.1. Also, Roche et al. (1999) obtained e ps 0. 

4. Comparison between model and data. Fit of ro 

and 7 

In the previous section, the only non-specified parameters 
are 7 and ro, on which the fluctuations strongly depend. 
These will be fitted in this section. Through the compari- 
son between model and data, we have a new method to ob- 
tain the two point correlation function, which does not re- 
quire explicitly the two-point angular correlation function, 
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Fig. 3. Log-log plot: fluctuations of the mean density of 
counts as a function of the radius of the circular regions. 
The error bar of the models is less than 1% . 
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Fig. 4. Modified x 2 -test. The contours show the values of 
7 and rj which are compatible with the data at 68% and 
95% C.L. 



do not suffer from edge problems, allows an estimation of 
the errors, and can be used even with incomplete surveys. 
Porciani & Giavalisco (2002) have also used galaxy counts 
for this purpose, but instead they derive the two-point an- 
gular correlation function and then insert it in the Limber 
equation (which assumes an exponential two point corre- 
lation function in all ranges), so it is halfway between our 
method and the standard one. 

Fig. shows the data and some model predictions. 
The numbers are in table The best fit is for ro = 6.66 
/i~ 1 Mpc, 7 ~ 1.79. Lower values of ro for the same 7, 
such as ro = 4.9 h~ 1 Mpc, give considerably less struc- 
ture than observed for the present sample of galaxies. For 
example, the two-point correlation function derived from 



6 Lopez-Corredoira & Betancort-Rijo: galaxy counts fluctuations 

Table 1. 8N/N as a function of 9q for the 2MASS-data and the models (same as Fig. |3J). 



0o (deg.) 


2MASS-data 


r = 6.66, 7 = 1.79 


r = 6.66, 7 = 1.79 Vr 


r = 4.5, 7 = 1.7 


r = 4.9, 7 = 1.8 


0.172 


0.869±0.004 


0.870 


0.873 


0.574 


0.678 


0.229 


0.785± 0.004 


0.777 


0.780 


0.521 


0.607 


0.286 


0.714± 0.004 


0.712 


0.715 


0.483 


0.557 


0.43 


0.599± 0.005 


0.606 


0.609 


0.422 


0.477 


0.573 


0.536± 0.006 


0.540 


0.544 


0.383 


0.427 


0.859 


0.461± 0.007 


0.457 


0.462 


0.336 


0.368 


1.146 


0.413± 0.009 


0.406 


0.411 


0.306 


0.330 


1.719 


0.353± 0.011 


0.342 


0.348 


0.270 


0.287 


2.292 


0.312± 0.013 


0.302 


0.309 


0.247 


0.259 


3.44 


0.254± 0.015 


0.253 


0.261 


0.217 


0.224 


4.58 


0.222± 0.018 


0.221 


0.231 


0.195 


0.200 


5.73 


0.197± 0.020 


0.198 


0.210 


0.178 


0.182 


6.88 


0.179± 0.022 


0.180 


0.194 


0.164 


0.167 


8.59 


0.156± 0.024 


0.158 


0.176 


0.146 


0.148 


11.46 


0.131± 0.027 


0.132 


0.154 


0.124 


0.125 


14.32 


0.106± 0.027 


0.114 


0.139 


0.107 


0.108 


17.19 


0.0895± 0.028 


0.0997 


0.127 


0.0946 


0.0955 


22.9 


0.0662± 0.029 


0.0798 


0.110 


0.0762 


0.0767 


28.6 


0.0535± 0.030 


0.0663 


0.0984 


0.0636 


0.0640 



0.21 
0.2 
0.19 
0.18 
0.17 
0.16 
< 0.15 
0.14 
0.13 
0.12 
0.11 
0.1 
0.09 



2 3 4 5 6 7 8 9 10 11 12 13 
# of resamplings 



Fig. 5. Value of A(Cy) as a function of the number of 
resamplings. 



APM optical survey (equal to that for the mentioned ha- 
los): 7 = 1.7, r = 4.5 fc^Mpc (Baugh 1996), shows a 
meaningful defect of fluctuations. This leads us to agree 
with Carlberg et al. (1997), that optically (blue) selected 
surveys appear to be significantly less correlated than K- 
selected galaxies. 

The values of the parameters 7 and r which are com- 
patible with the data are shown in Fig. 0] derived from 
a modified x 2 ~test for correlated data (Rubiho-Marti'n 
& Betancort-Rijo 2003, secc. 4) applied to the data for 
#0 < 7° (the other points have very large error bars and 
are not useful for constraining of the parameters; more- 



over, they are more dependent on the values of £ in the 
linear regime). 

In the simplest version of this test (Pi = Vi, see 
Rubino-Martm & Betancort-Rijo 2003), which for the 
present type of problem cannot be far from the best one 
(determined by an optimal set of Pi), the ordinary uncor- 
rected x 2 , Xord, and the number of degrees of freedom are 
rescaled by a certain factor, A, which is a function of the 
correlations: 



Xmod or a ' 

effective degrees of freedom 



(12) 



(13) 



where A depends on the correlations among data for all 
pairs of angular distances [see eq. (22) of Rubino-Martm 
& Betancort-Rijo 2003]. To obtain the confidence levels we 
proceed with Xmod' n mod like in an ordinary uncorrelated 
test with Xord> n ord The correlations among data, C\j = 
(jr(Qi)lT(Oj)} are calculated using resampling [jackknife 
resampling (Scranton et al. 2002, Mailer et al. 2003); 13 
resampling are generated and calculated in each one]. 
We find A(dj ) = 0.190 (in Fig. |S]we see how this value 
changes very little with the number of resamplings). 
The inferred values of the two parameters are: 



7 = 1.79 ± 0.02 



C.L.) [±0.07 (95% C.L.)], (14) 

(15) 



= 29.8 ± 0.3 (68% C.L.) [±1.1 (95% C.L.)] 



For 7 = 1.79, we derive r (z = 0) = 6.66 ± 0.04 (68% C. 
L.) h~ 1 Mpc (which is represented in Fig.[3J). 

These values are in excellent agreement with the es- 
timations of 7 in near- infrared surveys by other means: 
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the angular correlation function is at small scales propor- 
tional to 9~°- s in K-band (Baugh et al. 1996, Carlberg et 
al. 1997), which through Limber's equation gives 7 = 1.8. 
Kiimmel & Wagner (2000) give an angular correlation 
in K-band proportional to #-°- 98 ± - 15 which would give 
through Limber's equation 7 = 1.98 ± 0.15, again con- 
sistent with our value. Mailer et al. (2003) with 2MASS 
K-band data derives an angular correlation proportional 
to 6»-0-T6±o.04 which giveg 7 _ 17 q ± q Q4 5 V ery close to 

our preferred value. The value of 7 is also in agreement 
with some simulations of AC DM (Primack 2001, Fig. 1: 
7 = 1.6 — 1.7) and the exponents in optical surveys such 
as APM (7 = 1.7, Baugh 1996). 

We get an amplitude r^(z = 0) = 29.8 ± 0.3, more 
or less in agreement with the R-selected Las Campanas 
Redshift Survey which finds r^(z — 0.1) = 20.5 (Jing 
et al. 1998), equivalent [through eq. ijHJ with r comoving = 
?"phys.(l + z)] to Tq(z = 0) = 27.0. This is also in agreement 
with the extrapolation down to z — of the evolution of 
7-g in K-band selected galaxies at high redshift (see Fig. 8 
of Carlberg et al. 1997). However, the blue bands selected 
surveys give a lower amplitude of the correlations by a fac- 
tor of two (Carlberg et al. 1997, §6 and references therein). 
At high redshifts (z 3), this difference of amplitude be- 
tween K-band selected galaxies and blue-optical selected 
galaxies is even larger: a factor 3-4 (Daddi et al. 2003). 
This is consistent because the correlation length ro is not 
independent of the range of luminosities of the selected 
sample (Colombo & Bonometto 2001). Locally it follows 
roughly a dependence ro ~ 2.2 /i _1 Mpc +QADl with the 
scale D L = [< n > $(> L)]^ 1 / 3 (Colombo & Bonometto 
2001). 

Another important remark of Fig- El is that the fluctu- 
ations of galaxy counts may be accurately predicted by an 
expression © for the correlation. It must be noted that 
our method uses a frequentist statistical analysis that not 
merely determines the best correlation function but also 
tells us how good the best one is. The fact that the best 
values (of ro, 7) are within the 68% C.L. region tells us 
that expression JSJ) is as good as it can possibly be. 

No additional effect to large scale structure, for in- 
stance intergalactic extinction, is necessary to explain 
the variances. When Zwicky (1957, ch. 3) observed that 
&n(8o) I apoisson increased with 6*o even for large angles, 
he interpreted this result as a proof of the presence of an 
important intergalactic extinction. However, as we see in 
Fig. 01 the fact that 5N/N decreases slowly with increasing 
80 is due to the large scale structure itself, that produces 
extra fluctuations in large scales over those expected from 
a Poissonian distribution for scales of several degrees or 
tens of degrees. 

Summing up, we have tested that the ACDM model 
of large scale structure is consistent with the observed 
galaxy counts variance. This consistency is obtained for 
some particular values of the parameters of the two point 
correlation function, so a new method to obtain it indi- 
rectly without making use of the two-point angular cor- 
relation function. In the application of this method to 



2MASS K-band galaxies with wj- < 13.5 we get £(z 
II. r ■ 10 r'Mpc) = (29.8 iO.S)^ 17 ^ 02 (68% C.L.). 
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